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In this work, we introduce the hadronic loop contribution to explain the anomalous radiative transitions be- 
tween hi,{nP) and rjtimS), which was recently observed by the Belle Collaboration. Our calculation shows 
that the hadronic loop mechanism associated with these known decay mechanisms can explain why there exist 
anomalous radiative transitions between h),{nP) and ri),{mS). This study deepens our understanding of the decay 
mechanism of higher bottomonium radiative decays. 
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I. INTRODUCTION 

Besides reporting the evidence of a bottomonium 77^(25 ), 
the Belle Collaboration recently measured several branching 
ratios of the radiative transitions between htinP) (n = 1,2) 
and r]b{mS) (m = 1,2), i.e., !B[hb(lP) -^ rjb(lS)y] = (49.2 + 
5.7tli)%, S[hb{2P) -^ t]h{lS)y] = (22.3 + 3.8+^')% and 
!B[ht(2P) -^ 7]b{2S)y] = (47.5 ± 10.5!^;^)% [1]. As indicated 
in Ref. [1], there exists a discrepancy between these experi- 
mental values and theoretical expectations, where the branch- 
ing ratios [1] measured are a factor of 1.2 ~ 2.5 higher than 
the corresponding theoretical results given in Ref. [2]. This 
anomalous radiative transition between lib{nP) and rjbimS) 
stimulates us to propose a solution to this problem. 

In the past years, experimentalists have observed many 
anomalous hadronic decays and novel phenomena of higher 
charmonia and bottomonia, which include the excessive non- 
DD component of the inclusive t//{3nQ) [3, 4], the Okubo- 
Zweig-Iizuka (OZI) suppressed processes Xci -^ <^(j^, 4"P 
and double-OZI suppressed Xci — > ^0 [5], the Xcj radia- 
tive decays into a light vector meson [6, 7], the transitions of 
(/'(4040)/(/'(4160) into J/i/zi] [8-10], and anomalous large rates 
of e^e^ — > T(15,25)7r^:7r" near the peak of the T(55) res- 
onance [11]. Explaining these phenomena, some theoretical 
efforts have been rewarded [12-19], where the hadronic loop 
effects as an important QCD non-perturbative mechanism are 
introduced in these heavy quarkonium decays. Successfully 
explaining the observed experimental phenomena, we realize 
that the hadronic loop indeed plays a crucial role to the decays 
of heavy quarkonia. 

In addition, some discrepancies also appear between exper- 
imental measurements and theoretical predictions on the ra- 
diative transition between charmonia. The experimental mea- 
surements for J/i// -^ rjcj and if/' -^ flch'cJ [20] are much 
smaller than the nonrelativistic and relativistic Godrey-Isgur 
quark model predictions [2]. To alleviate the discrepancy 



between experimental measurements and naive quark model 
predictions, the meson loop contributions have been consid- 
ered in Refs. [21, 22] and good agreements with the experi- 
mental measurements have been archived. 

Along this line, in this work we adopt the hadronic loop 
effect to investigate the anomalous radiative transitions be- 
tween hb{nP) and rjbimS) observed by the Belle Collabora- 
tion [1]. In this study, we first need to answer whether the 
discrepancy between the experimental and theoretical results 
of these radiative decays can be alleviated by including the 
hadronic loop contribution added to the tree level diagram. If 
this answer is affirmative, the hadronic loop eff'ect, extensively 
applied to study intriguing higher charmonium and bottomo- 
nium decays recently observed, can be further tested, which 
will make our knowledge of non-perturbative QCD become 
abundant. In the next section, the details of the hadronic 
loop effects on hb{\P) —^ ?7i(15)y, hb(2P) -^ j]bilS)y, and 
hb{2P) —^ rih{2S)y will be explicitly presented. For the con- 
venience of our presentation, these transitions are abbreviated 
as hb{nP) -^ ribimS)y in the following sections. 

This paper is organized as follows. After the introduction, 
we illustrate the detailed formula of calculating the radiative 
transitions between hb{nP) and TjbimS). In Sec. Ill, the nu- 
merical results are given by comparing them with experimen- 
tal data. The last section is devoted to a short summary. 



II. RADIATIVE TRANSITIONS BETWEEN htinP) AND 
r]b{mS) 

In the naive quark model, the El transitions between P- 
wave and S -wave spin-singlets are written as [2] 



YQMihinP) -^ J]b{mS)y) = -ae^w^ \{'So \r\ 'Pi 



(1) 
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where a is the fine-structure constant and cq - -1/3 denotes 
the charge of the bottom quark in units of \e\. a> is the en- 
ergy carried by the emitted photon. The spatial matrix element 
{'5o \r\ ^P\) is relevant to the radial wave functions of initial 
and final heavy quarkonia. 

In Fig. 1, we give the typical quark-level diagrams depict- 
ing libinP) — > rjbimS )y in the naive quark model. In the naive 



quark model, a loop diagram is usually ignored when cal- 
culating the hi,{nP) — > rii,{mS)y transition. As indicated by 
Belle [1], the measured branching ratio of hh{nP) — > r]h(mS)y 
is not consistent with the result obtained in the naive quark 
model [2] which uses only Fig. 1 (a). Considering this situ- 
ation, we need to introduce the higher order diagram shown 
in Fig. 1 (b), where lii,{nP) first couples with two virtual bot- 
tom and anti-bottom mesons, which then turn into rii,{mS)y 
via exchanging bottom meson. However, it is rather difficult 
to calculate this kind of diagram Fig. 1 (b), so we would like 
to work in hadronic level diagrams. Here, we assume one to 
one correspondence between quark tree diagrams and meson 
tree diagrams as well as quark loop diagrams and meson loop 
diagrams. This is an assumption, but even now it is clear indi- 
rectly. It would be worthwhile to test it in future investigations 
when data will become more accurate. 




(a) 



FIG. 1: The typical quark-level diagrams for the hi,(nP) — > r]i,{mS)j 
transition. 



The hadronic tree level contribution to the El transition 
corresponding to Eq. (1) is given as follows. Assuming 
a point form factor of electromagnetic interaction and con- 
sidering the gauge invariance of the photon field, we can 
construct the Lorentz structure of the tree amplitude for the 
hb{nP) — > 77i(»j5')'y process as 



Mqm = gQMf^el^ [g^y - -2LJ!± 



Py ■ Phi, J 



(2) 



Here the effective coupling ^qm can be extracted by compar- 
ing the decay widths, Eqs. (2) and (1). The sign of ^qm is 
set to be positive in the present work. As indicated in the in- 
troduction, for hi,{nP) — > rib{mS)y transitions there exists a 
discrepancy between the experimental measurement by Belle 
[I] and the predictions by the naive quark model [2]. 

In the following, we introduce the meson loop contributions 
other than the tree contribution given by Eq. (2) when study- 
ing hi,{nP) — > rih{mS)y decays. Then we would like to answer 
the question whether the discrepancy can be alleviated by the 
hadronic loop effect or not. To calculate the diagrams shown 
in Fig. 2, the eff'ective Lagrangian approach is adopted, where 
we use the following Lagrangians constructed in heavy quark 
limit [23-27]: 






gh,B'Bh^,{B;B + b;b) + igh,B.B'^'' 

igni,B'BB'^(.d^JlhB - j]bd^,B) + H.C. 



dfihbvBlB'p, 



8nbB-B-s^ "/J 
ieAfjB- d B^ 



BlB^dfirjh, 
[ilsB'^ 



B-y 



^'"/"F^ySl^B- 



jgB'OBOy^""^F^vKliB°) + h.c] + ieA^ 



X [g-PBi- d B;'+g^''B:-d"B;^-g^"(fB:-B;' 



where S*g = daBl - dpB*^. In the heavy quark limit, the cou- 
plings among hi, j rib and the bottomed meson pair can be re- 
lated to two gauge coupling constants g\ and g2 by 



ghbB'B 
ghtB'B' 

gribB'B 
gilbB'B' 



-2gi ^Jmhi^niBinB' 

2g2^/m,,^mBmB', 
2g2mB' / y/ri^„ 



(3) 



where hb{nP) and rjbimS) are abbreviated as hi, and rjb, re- 
spectively. Thus, we need only to determine gi and g2 cor- 
responding to the couplings of hb{nP) and rjbimS), i.e., gi - 
- yl"^xw(nP)l'ilfxm(nP) and g2 = ^Jmry„s)/{2mBfr{mS))- The 
decay constants of P-wave bottomonium fxho(iP) and /^j„(2P) 
can be related to the radial wave function of bottomonia [28], 
i.e., /^,o(iP) = 265 MeV and /y„„(2/>) = 282 MeV. In ad- 
dition, /t(is) = 715 MeV and /t(25) = 482 MeV are de- 
cay constants of Y(15') and Y(25'), respectively, which are 
evaluated by their leptonic decay widths [29]. As for the 
_S*_Sy coupling constant, we adopt gfi.+B+y = 1.311 GeV"' and 
gB'OBy - -0.749 GeV"', both of which are estimated by the 
Ught- front quark model [30]. 
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FIG. 2: The hadronic triangle loop contributions to htinP) — > 
T}i,{mS)y. Replacing the charged bottom mesons with neutral bot- 
tom mesons in diagrams (c) and (d), we can obtain the remaining 
diagrams. 

Having these effective Lagrangians, we obtain the de- 
cay amplitudes corresponding to the diagrams listed in Fig. 
2 for the process hbipa) -^ 0*\pi)B'*\p2)]B'-*\q) -^ 



^b{P'i)y{pA), which read as 



gh,,B'Beh,,^, 



ignbB'Biipip + ipip) 



iee^(gai3(-iqv + ipiv) + gvp(-iqa) - gvai-ipi/i)) 



q - m- 



Mh = (/)■ 



r d^q 
J (2^ 



(4) 



gh,,B'B£h,,^. 



igntB-Biipip - iqp) 



ieeyv(-iq^ + /pj) 
1 



-gf'P + j/yjml, 1 



TO 9 9 

/'I - '«B- Pi - "^B 



q^ -rrp-g 



T{q\ 



(5) 



Mc = (/)■ 



J (27r: 






gh,,B'B^h„ 



- ginB'B'Se^pa/ii-ip'iXi/i) 



-^gB'BySp,iSTeyiiP4g'y - ipig^v)(-iqsg} 



+iq'gsi) 



2 2 

P2 - "4 






(6) 



/74 p 

T^ igh.B-B'Bpi^api-ipO^,, 



igni,B'Biip3A - iqi) 



9 J' 



TgB' By^etpST f y ( ip^g 



-iPlgy)i-iP2g^ + iPlgi) 



-g"'' + P"p\lml. 



-gi'^+i^Ai^ 



,f I, .a 
B' 



J 2 

q -m^ 



2 2 



T{q\ 



(7) 



/t4 r 

T^ igh,B'B-eppap{-iPQ)4„ 



■ gn„B'B'Se4,ST{-ip\){ip\) 



iee^(g,i£(-iqy + ip2v) 



+gv^{-iqA) - gvAi-ipi^)) 



-g'"'+p1p\lm'^. 



'1 "'B' 



9 2 2 '' 

P2 - nig. q"- - m-g. 



nq\ 



(8) 



where the form factor 'Fiq) is introduced to depict the internal 
structures and ofF-shell effects of the exchanged mesons as 
well as to remove the UV divergence in the loop integrals. In 
the present work, the form factor can be parameterized as [22] 



nq)^Y] 



(9) 



?; 



with qi - q,pi,p2, where the parameter A, can be 
parametrized as A, = m, + oAqcd with m, denoting the corre- 
sponding intermediate bottom meson mass and Aqcd = 220 



MeV. As illustrated in the caption of Fig. 2, there remain two 
more diagrams, whose amplitudes can be obtained by Aic and 
Aid by replacing the corresponding masses and the coupling 
constants. 




jfc(mA) -O 




'Ibi'nS) 



1 
B- B'- 

FIG. 3: Contact diagrams for hi, — » 7747. 

In order to keep the gauge invariance of photon fields intact, 
the contact diagrams given in Fig. 3 should also be included 
in our calculation. The effective Lagrangian describing vertex 
of TjhimS) interacting with B*B^*^y is 



-C, 



m,B'B<-'^y - gni,B'ByB*^ribA^tB + H.C. 

- ign„B-B-y>f'""'A^,B:BldiiTjb, 



where the coupling constants gri,,B'By and grj^B'B'y can be ob- 
tained by the corresponding ones of ■Cj^Bf-'iBi'h i-e., gq^B'By - 
eBginB'B and g^.B-By = eB-g,ji,B-B' with egieB') denoting the 
charge of B{B*) meson. By this effective Lagrangian, we 
obtain the amplitudes corresponding to the contact diagrams 
shown in Fig. 3, which are 



Mf = (/)' 



J (2^ 



ShB'B^,, 



q_ 

-gf^y + qftqv/m% 



gt]i,B'Byey 

1 



q'- - m-j^. 

j4 



iPo - q)- - nig 



rUq\ (10) 



Mo = {if 



I 7^-^ ighi,B'B'ep,,ai3(-ip^o)^i 



igtibB-B-ySvAetp^yiip^) 



-g"'^ + q"q''/ml. 



-gl'' + (pP,-qP){pl-q')lm 



q' 



TlM). (11) 



In the above amplitudes, another form factor Tconicf) should 
be introduced as a function of q^, which cannot be arbitrary in 
order to keep the gauge invariance. 

By summing over all decay amplitudes of the triangle 
hadronic loop diagrams and performing the loop integral, the 
general expression of amplitude for the triangle diagrams is 
parameterized as 



Wlxri = e^fft,, \gTng,iv - fln- 



PyvPhfi 



(12) 



Py ■ Phi, I 
On the other hand, the total amplitude of the contact diagram 



Mcon = Mf + Mg^ gCon44„8fiv 



(13) 



where the detailed deduction of the amphtudes for two contact 
diagrams is given in the Appendix. 



With the above preparation, we finally get the total de- 
cay amplitude of hhinP) — > rii,{mS)y considering only the 
hadronic loop eff'ect: 



MmL = Mfn + Wlco 



^7% 



igjn + gCon) gfiv - /xri 



PyvPhi,n 



Py ■ Ph, J 



r M V I PyvPl'b^|\ 

JMhe^£|,Ag^,v- - — — , 
\ Py- Ph„ I 



(14) 



which must be gauge invariant. In Eq. (14), we define /ml = 
gTii + gCon = /xri, where /ml or /m is obtained by calculating 
these triangle diagrams shown in Fig. 2, which is our main 
task in this work. 

As indicated above, in this work we will test whether in- 
cluding the hadronic loop contributions in hb{nP) — > rii,{mS)y 
can alleviate a discrepancy between experimental and theo- 
retical results. The total amplitudes for hb(nP) — > 7]h(mS)y 
transition is composed of AIqm given by Eq. (2) and A^ml, 
i.e.. 



Mot = MQM+e"''MMh 

I , i<t> f \ M V I PyvPhbH 

= [gQM + e ''/mlJ ety% \gtiv 

\ Py ' Ph 



(15) 



where (p is the phase angle between the amplitudes due to dif- 
ferent transition mechanisms. With these amplitudes, we can 
estimate the transition width as 



r{hb(nP) -> T]h{fnS)y) = Igqu + e'-^fi 



ML 



l2 ml - mj 
lAmnl 



(16) 



where niQ and m^ denote the masses of lib(nP) and rfhimS), 
respectively. 



III. NUMERICAL RESULTS 



TABLE L The y transition widths between hi,(nP){n = 1,2) and 
r]i,(mS){m = 1,2). The values of Tqm are estimated in the naive 
quark model [2]. Fml denotes the contributions from the meson loop, 
and their ranges are determined by parameter a. Lxot are coherent re- 
sults including the amplitudes of both quark model and meson loop 
contributions, which are dependent on parameter a and phase angle 
(f). The values of Lgxp are estimated from the corresponding branch- 
ing ratios from the Belle Collaboration [1] together with some theo- 
retical results; the details are given in the text. 

Initial Final Fqm [2] Fml Fxot Fsxp 
states states keV keV keV keV 
ht(lP) rib(iS) 37i) 0.1 -6.4 12.6-74.2 38.9-70.0 
ht(2P) rit(lS) 15.4 0.9-123.3 0-203.5 18.9-123.1 
ri,,{2S) 10.0 0.5-48.8 0-119.0 38.1-267.9 

Besides the coupUng constants mentioned in Sec. II, the 
masses involved in our calculation are ffJ/,,^(i/)) = 9899.1 
MeV, m,,„(2/>) = 10259.8 MeV, m;,,(is) = 9402.4 MeV, and 



f^niiiis) - 9999.0 MeV [1]. The properties of P and S wave 
singlet bottomonia are not well understood. At present, the 
full widths of hh{\P) and ht,(2P) are still not yet measured 
by experiment. Thus, we have to adopt some theoretical esti- 
mation, since we need to apply these resonance parameters to 
obtain the experimental partial widths of hh{nP) — » rib{mS)y, 
where the Belle Collaboration gave only the branching ra- 
tios of these decays [1]. Besides the radiative transition to 
lower bottomonia, hi,{\P) and hb{2P) dominantly annihilate 
into three gluons (ggg) or two gluons plus one photon (ygg). 
In Ref. [2], the widths of hb{lP) and hh{2P) decays into 
ggg and ygg were given, i.e., r{hb{lP) -> ggg) = 50.8 
keV, nhbilP) ^ ygg) = 1.6 keV, nhb{2P) ^ ggg) = 
50.2 keV, and Y(hb{2P) — > ggy) =1.6 keV. Thus we can 
obtain the branching ratio of hb{nP) — > rib{mS)y, which 
is S[hb{nP) -^ r]b(mS)y] = Th,(„p)^^,^(„s}yl'^l°l„py ^^^''^ 

^ht(nP) - 2m ^hi,(nP)^nh("^S)y + '^ h,,(nP)^ggg + '^h,,(nP)^ygg- With 

the partial decay widths of hh{nP) — > ggg /ygg from Ref. [2], 
we can roughly obtain the partial decay widths of libilP) -^ 
j]b(lS)y, hb(2P) ^ r]b{lS)y and lib(lP) ^ 7]h(2S)y as 
38.9 ~ 70.0 keV, 18.9 - 123.1 keV, and 38.1 - 267.9 keV, 
respectively. 

The lack of experimental measurement of the decay width 
of h„{\P) and hb(2P) leads to large uncertainties to the partial 
decay width of h bin P) — > yrjbimS). To compare the coherent 
partial decay width including the amplitudes of quark model 
and meson loop with those estimated from the branching ratio, 
we adopt a large parameter space, which is 1 ~ 4 for a and 
~ 27r for phase angle (p. The y transition widths between 
hb{nP) {« = 1,2) and rjbimS) {m =1,2) are present in Table I. 
The results from naive quark model are smaller than the lower 
limit of the experimental measurements, especially for hb{2P) 
radiative decay. The meson loop gives sizable contributions in 
the radiative decay of hb{lP)/hb(2P), and varies with parame- 
ter a. The discrepancies between experimental measurements 
and theoretical estimation in the naive quark model can be al- 
leviated. 

In Fig. 4, we show the decay widths of lib(lP) -^ r]b{lS)y, 
hb(2P) — > rib{lS)y, and hb(2P) — > rib(2S)y dependent on a 
in the definition of A, and a phase factor (p. The upper and 
lower limits of the experimental measurements are also pre- 
sented as solid curves. The parts sandwiched by the curves 
are the parameter spaces allowed by the experimental mea- 
surements. The comparison between our numerical and the 
experimental results indicates that the hadronic loop contribu- 
tions to hb{nP) -^ rib{mS)y indeed can alleviate the discrep- 
ancy as mentioned in Sec. I. Thus, anomalous radiative tran- 
sitions between hb(nP) and rjbimS) can be understood by in- 
troducing these hadronic loop diagrams. In addition, we also 
notice that there exist common a and (p ranges where the the- 
oretical values are consistent with the experimental data of 
hb{nP) — > T]h(mS)y (see Fig. 4 for more details). This phe- 
nomenon also reflects the similarity existing in the lib{lP) — > 
T]b{lS)y, hb{2P) -> r]bilS)y, and libi2P) -> T]b{2S)y decays. 

At present, there exists a large experimental range for these 
discussed hb{nP) — > rih{mS)y transitions [1]. We expect 
that future experiments can give a more precise measurement, 
which will be useful to further constrain our parameter range. 




^h{2P)^r,t{2Sh (keV) 




I 110 
' 100 
90.0 
80.0 
70.0 
60.0 
50.0 
40.0 
30.0 
I 20.0 
10.0 
0.00 



FIG. 4: (color online). The a and dependence of the partial de- 
cay widths ofhtilP) -^ r]b(lS)j, h,,{2P) -^ ri,,{lS)y, and h{2P) -^ 
Tib(2S)-y con'esponding to the diagrams from top to bottom, respec- 
tively. Here, the black curves are experimental lower and upper limits 
of the corresponding decay width given by Belle [1]. Units of values 
for decay widths are in keV. 



IV. SHORT SUMMARY 

Being stimulated by Belle's observation of anomalous ra- 
diative transitions between hi,{nP) and Tii,{mS), in this work 
we study hi,{nP) — > rii,{mS )y transitions by introducing the 
hadronic loop contributions. As shown in our numerical re- 
sults, the hadronic loop contributions play an important role 
to get consistent results with the experimental data [1]. These 
phenomena also show that the hadronic loop effects can be 
a universal mechanism existing in the charmonium and bott- 
tomonium hadronic and radiative decays, since there have 
been some theoretical efforts [12-19] in this direction. 

With more and more accumulation of experimental data, 
some novel phenomena of the decays of higher charmonium 
and bottomonium have been revealed. If inclusion of the 
hadronic loops is a universal non-perturbative QCD effect, we 
believe that the hadronic loop mechanism should be further 
tested in the future by comparing the results with more ex- 
perimental observations, which is an intriguing and fruitful 
research field. 
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Appendix A: Contributions from Contact diagrams 

As for the amplitude Aif, one has 



Mf - ^el 



J (2;r)n.2_, 



Tconiq') 



[q~ - ml.] [ipo - q)^ - «4] 
" 4^k {foSiiv + /i/^opPOv) = fo^elgi^y. (Al) 

For the ampUtude Aig, it is in the form 



Mg oc epf,a/3i:>Q^ SvAe<i,e- 






d^q -g"^+q"q^/ml. 



ipu - qf - ml. 



^Co„(?') 



X 1 /{[q^ - ml. ] [{po - qf - m^T^.M^) 

= ep^aplfAe,Ae^e;pl[gdggr''^ + gdgPoPor''" 

+ gilpoPoPoPoT'"') 

where [gg]"^'''^ indicates the terms in which these four Lorentz 
indices are shared with two metric tensors, and [gpopo]'"'^^'^ 



and [popopopoT'^^'^ are defined in the same way. It is easy to 
find that the term proportional to g2 vanishes after contracting 
all the Lorentz indices. The symbol [gpopo]"''^'^ includes two 
different cases. One is that a and/? are the Lorentz indices of 
the metric tensor, i.e. g"^pf^plf. The other is that at least one 
of a and /? is the Lorentz index of momentum po; i.e., these 
terms are proportional to p"^, pl^ or PqPq. After contracting all 
the Lorentz indices, one can find that all these terms in two 
cases result to zero. Then for Mg only the terms proportional 
to gQ survive, that is, 

Mg oc Epf,appf^€^f^eyAe4:e:ypt,[ggf'^^'' 
°^ fSpt^apf^^h.^vAe^eyPlf'^g^" 
= 2po ■ Pi^t^/yg^tv - 2p3 ■ ei,,po ■ £y. (A2) 

In the initial state rest frame, pg = mo and /?o = 0, but for 
a real photon, the polarization vector has efj = in both the 
Coulomb and axial gauge, and , thus, we have po ■ ey -Q, and 
the amplitude Mg should be in the form 



Then according to Eqs. (Al) and (A3), we have 



(A3) 



A^Con = Mf + Mg = gCon^/ygt^y (^4) 
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